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' Abstract. We present a new determination of the local temperature function of X-ray 
' clusters using a sample of X-ray clusters with fluxes above 2.2 10~^^erg/s/cm-^ in the 
[0.1 — 2.4] keV band, most of these clusters come from the Abcll XBAC's sample to which 
■ a handfuU of known non-Abell clusters has been added. We estimate this sample to be 
85% complete, and should therefore provide an useful estimation of the present-day num- 
^ , ber density of clusters. Comprising fifty clusters for which the temperature information 
^ ,: is available, it is the largest complete sample of this kind. It is therefore expected to 
significantly improve the estimation of the temperature distribution function of clusters. 



We find that the resulting temperature function is higher than previous estimations, but 
it agrees with the temperature distribution fimction inferred from the BCS and RASS 
luminosity function (Ebeling et al., 1997; Do Grandi et al. 1999a). We have used this 
sample to constrain the amplitude of the matter fluctuations on cluster's scale of 
8v'T1o~^ft,^^Mpc, assuming a mass-temperature relation based on recent numerical sim- 
ulations. We find (Tc ^ 0.6 ± 0.02 for an fio = 1 model (for which Uc = cts). Our sample 
provides an useful reference at z ~ to use in the application of the cosmological test 
based on the evolution of X-ray clusters abundance (Oukbir & Blanchard 1992, 1997). 
We have therefore estimated the temperature distribution function at z — 0.33 using 
Henry's sample of high-z X-ray clusters (Henry, 1997; hereafter H97) and performed a 
preliminary estimate of flo- We find that the abundance of clusters at z — 0.33 is signif- 
icantly smaller, by a factor larger than 2, which shows that the EMSS sample provides 
strong evidence for evolution of the cluster abundance. A likelihood analysis leads to a 
rather high value of the mean density parameter of the universe: f^o — 0.92^q'215 (Icr) 
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for open universes and Oq = 0-865lo!245 universes, which is consistent with a 

previous independent estimation based on the full EMSS sample by Sadat et al.(1998). 
Some systematic uncertainties which could alter this result are briefly discussed. 

Key words: Cosmology: observations - Cosmology: large-scale structure of the Universe 

- Galaxic's: clustors: gcnc^ral 



1. Introduction 

Clusters arc believed to be the largest virialized concentrations of dark matter. Therefore, 
they offer privileged regions for studying dark matter distribution on large scales in the 
universe. X-ray and weak lensing mass measurements have been added to the traditional 
mass estimates based on optical velocity dispersions. Weak lensing analyses are still in 
their infancy and at present there exists no sample of weak lensing observations large 
enough to establish a mass function. The velocity distribution function of clusters could 
be used to derive the mass function. However, as Evrard pointed out (Evrard, f989) 
the error on individual measurements can introduce a significant overestimate. Further- 
more, velocity measurements at least in the case of distant clusters, can be corrupted 
by projection effects that might be difficult to remove in practice (Frenk et al., 1990). 
Moreover, Sadat et al. (1998) have shown that X-ray temperatures of some of the CNOC 
clusters show a significant difference with what is expected from their velocity dispersion 
measurements. For these reasons, it has been argucid that the X-ray temperature is a bet- 
ter indicator of cluster mass. Numerical simulations have greatly helped to understand 
the relation between X ray temperature and the mass, and useful constraints have been 
placed on the amplitude and the shape of the spectrum of mass density fluctuations. Still, 
the size of the samples of X-ray clusters homogeneously selected is limited: 25 clusters in 
the Henry and Arnaud (1991, hereafter HA91) sample and 30 clusters in Markevitch's 
(1998). ROSAT selected clusters samples have significantly improved the situation in 
this domain. In order to achieve stronger constraints on theoretical models, it will be 
necessary to obtain more temperature measurements, a vast program that will probably 
be possible with the next generation of X-ray satellites such as AXAF and XMM. 

The cluster X-ray temperature function is a powerful tool for cosmology. Provided 
that the mass-temperature relation is reasonably well known, the Press and Schechter 
(1974) formalism allows one to constrain the amplitude and the shape of the power spec- 
trum for a given cosmologieal background density (see Bartlett, 1997, for a recent review 
on the subject). Since there is a nearly complete degeneracy between the amplitude 
of the fluctuations and the mean density of the imiverse, the cosmologieal parameter 
r^o cannot be determined solely from the local temperature function. However, the 
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evolution of this temperature distribution function, once normalized to the present day 
cluster abundance, varies significantly with offering an interesting new cosmological 
probe (Oukbir & Blanchard 1992, hereafter OB92; Hattori & Matsuzawa, 1995). This 
test has received considerable attention in recent years (Donahue, 1996; Carlberg et al., 
1997; H97; Oukbir & Blanchard 1997, OB97, hereafter; Blanchard & Bartlett, 1998; Eke 
et al., 1996, 1998; Sadat et al., 1998; Viana & Liddle, 1996, 1999a; Blanchard et al., 
1999; Donahue et al., 1999; Donahue and Voit, 1999; Reichart et al., 1999). Variants 
have been proposed using the Sunyaev-Zeldovich (Barbosa et al, 1996) and weak lensing 
effects (Kruse & Schneider, 1999). Modeling the redshift distribution of the EMSS X-ray 
selected sample (OB97) given the absence of any significant evolution in the — 
relation (Mushotsky & Sharf, 1997; Sadat et al., 1998), seems to favor high value of the 
density parameter (Sadat et al, 1998; Reichart et al, 1999). Modeling the RDCS redshift 
distribution leads to consistent results (Borgani et al, 1998). A more direct estimate, 
free from any consideration on the possible evolution in the — relation, could 
be obtained from the measurement of the evolution of the temperature distribution 
function, which in turn requires a good knowledge of the selection function of the 
sample of clusters. The first sample of X-ray selected clusters at non-zero redshift with 
measured temperatures has recently become available (H97) and has led to an apparent 
median value of fig in the range 0.2 — 0.6 (H97; Eke et al, 1998), although higher values 
were found by Viana and Liddle (1999a) and Blanchard et al. (1999). It has been argued 
that current data are not good enough to allow a reliable estimate of the mean density 
of the universe from such techniques (Colafrancesco et al., 1997), a conclusion that 
appears wise given that the local sample of X-ray clusters used up to now is that of 
HA91 which contains only 25 clusters while the high-redshift sample comprises only 10 
clusters with moderately accurate temperature measurements (recently, the redshift of 
one cluster in the sample has been revised, reducing the sample to 9 clusters, Donahue 
et al., 1999). However, the fact that some conclusions have already been drawn, even 
if too optimistic regarding possible systematics, demonstrates the power of this test: 
clearly ten more clusters or so at high redshift and an accurate determination of the 
local temperature distribution function would provide a very robust determination of the 
mean density of the universe. In fact, going to high redshift makes a dramatic difference 
in the abimdancc of hot clusters (OB92), and the existence of few EMSS clusters with a 
high temperature has been argued to already provide a strong evidence for a low den- 
sity universe (Donahue 1996; Bahcall & Fan, 1998; Donahue et al., 1998; Eke et al., 1998). 

Until recently, the X-ray cluster temperature distribution function was inferred 
from catalogs built from the HEAOl survey in the 2-10 keV band. The need for a new 
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estimation of the temperature distribution function has been recognized. Markevitch 
(1998, hereafter M98) has provided such a new estimate based on a sample of X-ray 
clusters for which ROSAT fluxes were available. Most of his clusters come from the 
XBACS sample of bright Abell clusters (Ebeling et al., 1996). Temperatures are derived 
from ASCA observations, and for this reason clusters at redshift smaller than 0.04 were 
not considered. This means that the sample is not a strictly flux limited sample. M98 
corrected both the fluxes and temperatures of his sample for the presence of cooling 
flows in the central regions. Cooling flows are believed to be an important source of 
dispersion in the — relation (Fabian et al., 1994; M98; Arnaud & Evrard, 1999), 
but have negligible effect on the temperature distribution function. 

In the following, we present the analysis of a new sample containing flfty clusters 
essentially based on the XBACS sample (Ebeling et al., 1996), using a flux limit sim- 
ilar to that of M98. Note that M98 considered his sample to be reasonably complete. 
In the present analysis, the effect of temperature error measurements is explicitly taken 
into account using a Bayesian correction. Our sample does not require any correction 
for redshift incompleteness. We also do not correct for cooling flows in order to allow a 
direct comparison between our sample and those at high redshift. The resulting temper- 
ature distribution function we obtain is smooth and can be used for useful comparison 
with theoretical models. We have performed such a comparison in order to constrain the 
density fluctuation power spectrum with different values of the cosmological background 
density Oq- This paper is organized as follows: section 2 summarizes the results from pre- 
vious estimates of the local temperature function. In section 3 we give a brief description 
of our sample. In section 4, we present our method to calculate the local temperature 
function as well as a determination of the distribution function of the cluster abundance 
estimator. We also discuss how our temperature function compares with previous results. 
Section 5 presents the flrst application of this new temperature distribution function to 
constrain cosmological parameters. 

2. Previous estimates 

In order to determine the temperature funetion, it is important to have a well controlled 
sample of clusters. The knowledge of the selection function is therefore critical. For in- 
stance, the completeness of the Edge et al. (1990) sample (hereafter E90) was questioned 
by HA91, and in fact the difference between the E90 and HA91 temperature fmictions 
is striking. We therefore re-examine these two samples, although the revised version of 
HA91 appears to be closer to E90. Rather than plotting the temperature fmiction, we 
directly examine the histogram of the number of clusters in temperature bins for each 
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Fig. 1. The distribution of clusters with revised temperatures in the Henry & Arnaud 
sample, HA91 (left panel). The right panel shows the histogram of the 21 clusters in the 
Edge et al. (1990) sample, not included in HA91. 



sample. In Figure we show the histogram of the revised temperatures in the HA91 
sample, while in Figure |l|b, we give the corresponding histogram for the E90 clusters not 
present in the HA91 sample; this latter will be referred to as the E-HA sample. Com- 
parison of the two histograms clearly reveals a striking difference: there is an apparent 
deficit of clusters with T ~ 4 — 6 keV in the HA sample as compared to E-HA. Obvi- 
ously, the smaller number of 4 keV clusters present in the HA sample leads to a lower 
number density. This is of critical importance in the theoretical modeling, because the 
abundance of 4 keV clusters is crucial in determining the amplitude of the fluctuations 
on a 8/i^^Mpc scale (for Hq = 1). The E-HA set is astonishingly different: the number 
of clusters reaches a maximum at T = 4 keV, and then tends to decrease monotonically 
for higher T, meaning that there are relatively fewer hot clusters in the E-HA data set 
(4 clusters from the E90 sample with poor temperature measurements are not taken into 
account in E-HA, which adds 3 clusters to the region T > 8 keV, but this does not 
change the shape of the histogram) . Although statistical analysis of such a small sample 
is a hazardous exercise, this difference is highly significant: assuming that the number 
of clusters found in the HA sample between 4 and 6 keV (4 clusters), the probability to 
have ten clusters or more in the E-HA sample in the same range is less than lO^'^. When 
dealing with small samples, one must of course be careful of spurious noise introduced by 
Poisson statistics. However, using different tests we always find this difference significant 
at a confidence level greater than 95%. It is difficult to understand the origin of such a 
difference: both samples are X-ray selected and incompleteness is expected to be nearly 
independent of temperature. 

It is possible that the existence of large-scale fluctuations, as can be seen from the 
visual appearance of the clusters distribution or from the long-range correlation function, 
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indicates that noise due to a finite number of objects may be larger than that expected 
from Poisson noise. 

The comparison we have presented in Fig. |l| clearly calls for a larger, complete flux 
limited sample strictly X-ray selected. However, only the Rosat All Sky Survey (RASS) 
would allow the construction of such a sample. The Bright Cluster Survey BCS (Ebeling 
et al., 1998) and the RASS South (Dc Grandi et al., 1999b) have just become available at 
the time of writing, but the lack of temperature information for a significant fraction of 
these clusters prevents us from using these samples to usefully estimate the temperature 
distribution function. Recently, M98 provided a new local temperature function based on 
the XBACS sample (Ebeling et al, 1996). XBACS is an essentially complete, flux- limited 
X-ray sample of Abell clusters. In principle, the restriction to Abell clusters implies that 
the sample is possibly incomplete. There are three possible sources of incompleteness: 
i) low mass (i.e., low temperature) clusters may be missing because they are not rich 
enough to be selected by the Abell criteria; ii) distant clusters could be missed because 
at large redshifts clusters hardly meet the Abell criteria; iii) there is a significant difference 
between the optical based Abell criteria and those based on X-ray fiuxes. M98 restricted 
his sample to the redshift range 0.04 < z < 0.09, the lower limit being imposed by 
ASCA. This redshift restriction eliminates intrinsically faint clusters from his sample. 
M98 corrected for this incompleteness by using a weight to compensate for the selection 
function. The resulting temperature function is noticeably higher than previous estimates, 
although it is claimed to be consistent with them. However, the correction procedure 
could introduce a systematic error in the derived abundance of cool clusters as we will 
see in the next section. 

3. The present sample 

As it is clearly an advantage to have a survey which is not redshift restricted, we have 
collected existing information on clusters with fx > 2.210^^^ crg/s/cm^ in the ROSAT 
band 0.1-2.4 keV (and with > 20°). In practice, most of the clusters were obtained 
from the XBACS sample complemented by a handful of non- Abell clusters which satisfy 
our fiux requirement (MKW3s, Virgo, C10336 and 0422-09). For all these clusters, tem- 
perature measurements were found in the literature. In the following analysis, the redshift 
range is still restricted to z < 0.1, because of the possible substantial incompleteness of 
the sample at higher redshifts. Our final sample contains fifty clusters. It is important to 
examine its completeness. 

A classical test for checking the homogeneity of a survey is to use the V/Vm test 
(Schmidt, 1968). Applying this test, we find a mean value of V/Vm = 0.490, consistent 
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with a homogenous sample. The plot of V/Vm versus temperature given in Figure ^ 
does not reveal any unexpected trends. The recent availability of the BCS and RASS 
samples allows us to check this in a more quantitative way: at our chosen flux limit, 
four clusters in the BCS and RASSl are not present in our sample (RXJ0419, IIZwlOS, 
Z5029, NGC1550). None of these (non-Abell) clusters has a measured temperature, but 
from their luminosity, only two of them would have a temperature greater than 4 keV. 
This is consistent with our sample being complete at a level of 95% for T > 4 keV, while 
the completeness of the BCS and RASS is estimated to be better than 90%. Therefore, 
our estimated completeness is 85% and we will treat it as an essentially complete, X- 
ray selected ensemble of clusters, at least for clusters with temperature greater than 4 
keV. It is the largest homogenous sample of X-ray selected clusters currently available 
for the estimation of the temperature distribution function. Notice that the possible 
incompleteness of our sample would imply that the inferred temperature function could 
be an under- estimate. 

At redshifts z > 0.04, our sample is similar to M98. The temperature histograms of 
both M98 sample and ours are shown in Figure ||. As one can see, the two histograms 
are nearly identical on the hot end, but differ on the cool end (for T < 5 keV). This 
difference is not surprising given the restricted redshift range investigated by M98. 
However, it is important to investigate the temperature function on the cool end since 
the number of clusters in the two samples differs noticeably. Although M98 corrected 
for the incompleteness of his sample, using a larger sample will result in a statistically 
better estimate with reduced sensitivity to systematics. Such a sample has two main 
advantages: with a selection in the ROSAT band, we are much less sensitive to the 
selection function of the 2.-10. keV band, which obviously favors hot clusters (even if 
in principle the effect of the selection function can be corrected for), and the number 
of clusters in our sample is twice as large as the HA91 sample, leading to better 
statistics especially on the cool end (T ^ 4 keV). We have examined the temperature 
histograms for clusters with fluxes in the range 2.2 — 4.0 10^^^ erg/s/cm^ and those 
with fx > 4.0 10^^^ erg/s/cm^. We flnd that these two histograms differ noticeably 
in the same sense as the HA91 and E-HA samples as discussed previously. The reason 
of this difference is unclear, but it is clearly not due to the 2.-10 keV selection only. 
Inspection of the redshift distribution does not provide evidence that this could be 
due to large-scale structures (although the BCS reveals significant fluctuations around 
z - 0.05). 



4. The Determination of the Temperature Function 



Blanchard, Sadat, Bartlett & Le Dour: Cluster Temperature Function 



1.0 



0.8 - 



0.6 - 



E 
> 



0.4 - 



0.2 - 



0.0 



o o 

^ 



o o 
o 



6 

T (keV) 



10 



Fig. 2. V/Vm distribution for clusters in our sample as a function of temperature. No 
sign of incompleteness can be seen from this plot for clusters with T > 3.5kcV. 





a) b) 

Fig. 3. Same as in Fig. ^ for clusters selected in the ROSAT band. Figure a) represents 
the distribution of clusters in Markevitch's sample, b) represents the histogram of the 
50 clusters in our sample (XBACS clusters with flux greater than 2.2 lO^^^erg/s/cm^ 
and redshift z < 0.1). The shapes clearly differ for temperature smaller than 5 keV, 
while they are nearly identical above. This difference is mainly due to the fact that the 
Markevitch sample is restricted to redshift > 0.04. 



4-.1. The method 

Estimation of the luminosity distribution function (f>{L) from a sample for which the 
selection function is well known is relatively straightforward: the observed number of 
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clusters in the sample Af in the luminosity range {L ~ L + AL) (with AL <C L) is the 
realization of a stochastic Poisson process with a mean value (j){L)V{L)AL, where V{L) 
is the volume searched for objects with luminosity L. When the selection function s{z) 
is that of a strictly flux limited sample, V{L) is just the volume Vm corresponding to 
the maximum distance at which the cluster would have been detected. If the selection 
function s{z) is more complex, then this volume can always be computed as an integral: 

V{L) = / s{z)dV{z) 







An unbiased estimation of 0(L)AL is therefore given by J\f/V{L). In the following, 
when B is an estimator of some quantity A, we use the notation below: 

A = B 

According to this notation: 
cj){L)AL = N/V{L) = ^/V 

where the summation is performed over all the relevant clusters in the sample. This 
immediately provides an estimator of the integrated luminosity function: 

0(> L)^J2 

Li>L 

We can now recover the standard estimator of the temperature distribution function: 

r + co 

N(T)AT^ / (j){L,T)ATdL = y2'^/ViT^) (1) 
Jo 

where the summation is performed over all clusters in the sample in the relevant temper- 
ature range. It is important to notice that this estimator is unbiased and fully accounts 
for the intrinsic dispersion in the temperature-luminosity relation. This procedure may 
not always be used in practice, if for some region of the parameter space in the sample 
V{L) = 0. This is the case for instance, in a flux limited sample which is further re- 
stricted to a redshift range. Because of this further restriction, clusters fainter than some 
intrinsic luminosity may simply not be present in the sample, consequently leading to 
V{L) = 0. In this case, the estimation from Eq. (|^) is clearly not appropriate anymore. 
The contribution of the missing clusters can however be estimated, provided the prob- 
ability distribution of clusters at a given temperature p{L\T) is known. The observed 
number of clusters is still a Poisson realization of a process whose mean is given by: 

r + co r+oo 

/ (l){L,T)ATV{L)dL^ I N{T)ATp{L\T)V{L)dL 

Jo JQ 



N{T) / p{L\T)ATV{L)dL = N{T)ATVs{T) 

JQ 



(2) 
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The last equality defined Vs{T), the effective average volume per cluster of given tem- 
perature in the survey. Therefore: 

iV(T)AT = ^l/y,(ri) 

i 

This estimator has been used by several authors (Eke et al., 1998; M98). It is unbiased 
provided that Vs (T) is known exactly. In practice, this is of course not the case, since the 
statistical description of the temperature-luminosity relation and the distribution p(L|T) 
are not perfectly known, which leads to a source of systematic uncertainty. Consequently, 
a specific choice of p(L|T) implies a bias in the above estimator. Furthermore, as fewer 
clusters are actually present in the sample, the Poisson noise increases. It is clear there- 
fore, that this approach should be avoided when possible, unless one can be confident 
that the resulting uncertainties remain small. 

A second important issue concerns errors on temperature measurements: when the 
true underlying cluster distribution is very steep, positive errors will move many more 
clusters upwards in temperature than negative errors move hotter clusters downwards, 
just because of the difference in the intrinsic abundances. This effect has been pointed out 
by Evrard (1989) for velocity dispersion measurements. A mean statistical correction was 
applied in the modeling by Eke et al.(1998), while Viana and Liddle (1999a) estimated the 
magnitude of this effect using a bootstrap resampling technique on Henry's sample. It is 
important to realize that when the errors on individual measurements are not identical, a 
mean correction is not necessarily sufficient, because the weight 1/Vi varies significantly 
from one cluster to another. A further problem is that, in principle, the error on the 
temperature might correlate with the apparent flux of a cluster, as fainter clusters will 
have lower signal-to-noise and will be those having the larger 1/Vj as well. In such cases, 
only a Monte-Carlo reproducing the exact conditions of the observations including the 
different integration times for the different sources, would in principle allow a complete 
separation of the various effects of temperature measurement errors. 

For the local temperature distribution function, the errors are small enough that 
the above issue is not a major problem. However, due to larger temperature errors this 
might be a more critical point for the high-redshift sample. The bootstrap method used 
by Viana and Liddle (1999a) is certainly well adapted. Another possible way to solve this 
problem is to use a Bayesian approach: we can take as a prior the distribution function of 
X-ray temperatures assuming N{T) oc T~", with a w 5. Given that a cluster is observed 
with a measured temperature Tq and assuming that the errors are log- normal distributed 
with a dispersion a (corresponding the observed la uncertainty), the a posteriori prob- 
ability that the actual temperature of the cluster be T is given by: 

p{T\To) cx N{T) exp(-(ln(r) - ln{To)f /2.a^) (3) 
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this distribution is also log-normal with the same dispersion but with the most likely 
temperature T shifted compared to Tq: 

T = roexp(-acr^) (4) 

The above formula may appear strange, because it seems to suggest that a systematic 
bias exists in the temperature measurements. This is actually not the case, since the 
observed number of clusters within the sample is rather constant among the different 
temperature bins. Only the estimated density is biased. 

Another effect to take into account for cosmologieal applications (see section 5.3) 
when one wishes to relate the mass function to the temperature distribution function, is 
the dispersion in the M — T relation. The importance of this effect has been discussed by 
Eke et al (1998). They assumed an intrinsic dispersion of 20%, but neglected the effect 
of the errors on the actual measurements considering them to be smaller. However, as 
the effect goes as the square of the error, one should be cautious when dealing with these 
corrections; in practice, a 20% error (or dispersion in the M — T relation) begins to make 
a significant difference in the inferred abundance of clusters, while an amplitude of 10% 
leads to a change that is essentially negligible. In their numerical simulations, Bryan and 
Norman (1998a) found an intrinsic dispersion of the order of 10% in the M — T relation, 
twice smaller than the dispersion adopted by Eke et al.(1998) (w 20%). In order to 
understand the bias introduced by these effects (errors on temperature measurements 
and dispersion in the M — T relation), the Bayesian correction is illuminating: as we have 
seen, such a correction is equivalent to a modification of the observed temperature, that 
is to say, the mass-temperature relation in the modeling (this is only an approximation 
because the Bayesian term depends on the actual slope of the N(T), which might vary 
with T and redshift, but this is a second order effect). The corrections to implement 
for observational errors arc slightly different in nature though: errors vary from one 
cluster to another, and also may vary in a systematic way, for instance accordingly to 
the apparent flux. Furthermore, since the correcting factor varies in a nonlinear fashion, 
a mean correction is inadequate. Finally, it is worth noticing that errors for high redshift 
clusters are larger than for the low redshift clusters, and therefore a larger correction is 
necessary. The practical implication of such a change in the modelling will be discussed 
is section 5. 
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4-. 2. Temperature Distribution Function 

We compute the integrated temperature distribution from our sample following the 
method described in the previous section, i.e. using the estimator: 

N{> T)^ J2 

Ti>T 

The resulting temperature distribution function is plotted as the continuous thin line in 
Figure ^. A smoothed version is also given. We have also checked that using temperature 
when the central emission is removed makes negligible difference. For comparison, we 
have also estimated the temperature distribution function from the HA91 sample, using 
available updated temperature measurements and the standard estimator (Eq. |l|) without 
the Bayesian correction factor. This temperature distribution is shown as the dotted line 
m Fig. I As one can see, the difference is quite noticeable, in particular around 4 keV, 
where it is roughly a factor of two. We also show the temperature distribution inferred 
from the BCS luminosity function (Ebeling et al, 1997) assuming the ~ relation 
as given by OB97. A convenient fit of our temperature distribution function is: 

7V(> T) = 1.4210~^/i^Mpc"^(r/4.keV)"^-^exp(-(r-4.keV)/2.6) (5) 

which is plotted in Figure |[ We also show the fit to the temperature distribution 
function derived by M98. As one can see, we find a larger abundance of cool clusters 
(r < 4 keV) than M98 while we have a nearly perfect agreement with the temperature 
distribution function inferred from the BCS luminosity function. 

The difference between our temperature distribution function and previous estima- 
tions is striking. The difference with M98 can be easily understood - his sample being 
restricted to the redshift range z > 0.04, a significant fraction of faint (and cool) clusters 
are missing. He has corrected for this incompleteness by a weighting scheme based on 
an estimated dispersion in the L^. — relation and following the method discussed in 
section 4.1. However, in order to properly evaluate such a correction, one needs an accu- 
rate estimate of the dispersion in the Lx — Tx relation, which can only be obtained from 
a flux-limited sample. The difference with the original HA91 sample is more difficult to 
understand. One may argue that the difference comes essentially from the band selection 
(2 — 10 keV ), however as we have already mentioned in section 2, comparison with the 
E-HA sample suggests another origin. We have therefore divided our sample into two 
equal sub-samples, corresponding to the brightest and the faintest clusters (in apparent 
flux). Again, since both samples are X-ray selected (in the same band), they should be 
statistically equivalent. The corresponding temperature distribution functions are shown 
in Figure ^. The bright sub-sample leads to a temperature distribution function close to 
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Fig. 4. The integrated temperature distribution function inferred from our sample is 
given by the continuous (red) thin hne. A smoothed version is also given (thick line). 
The dotted line is the same quantity for the original HA91 sample estimated by means 
of Eq. |l]. The dashed (green) line is derived from the BCS luminosity function. The error 
bars represents the 68% interval for the distribution of the estimator from the bootstrap. 

the one based on the HA91 sample. This is not surprising on the high temperature end, 
since the clusters are almost the same in these two samples. The fact that the tempera- 
ture distribution functions are almost the same for the cooler clusters as well, indicates 
that the original HA91 sample does not suffer from any significant bias and that the 
effects of the chosen band are relatively small. On the other hand, the difference between 
the faint and the bright samples is quite significant, and the relative abundances of hot 
and cool clusters is noticeably different. As the difference is marginally 1 cr, these two 
samples could be considered as representing two realizations which are slightly in the 
tail of the distribution, without being really suspicious. Nevertheless, such a difference 
have some implications in practice as we will see in the next section. The fact that our 
temperature distribution function is in good agreement with that inferred from the BCS 
luminosity function is rather encouraging and confirms the absence of any systematics 
that would undermine our analysis. 

Differences such those we have found have significant consequences when constraining 
the spectrum of the fluctuations from the temperature distribution function: a sample 
in which the abundance of cool clusters is underestimated leads to a too flat spectrum 
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Fig. 5. A fit to our temperature distribution function is given by the continuous (red) 
tliick line. The two thick dashed lines correspond to two sub-samples: the short dashed 
line corresponds to the bright one (2.210""'^^ > > 4.10"^"'^ erg/s/cm^), the long dashed 
corresponds to the faint sub-sample {sx > 2.210^^^ erg/s/cm^). The long-dashed-3- 
dotted (green) line is the temperature distribution function from Markevitch (1998). 

and the amplitude of fluctuations at cluster scales will be also underestimated. All these 
effects enter as a source of systematic uncertainty in the estimation of Qq. 



4-3. Distribution function of the estimator and Confidence intervals 

It is useful to have a statistical description as comprehensive as possible in order to 
derive strong constraints on models. The number M of clusters in a given temperature 
range can be regarded as a Poisson realization of mean /q °° (j){L) V(L) A L. However, 
the understanding of this quantity requires good information on the — relation: its 
shape, its normalization and the distribution function around it. In order to determine the 
expected distribution of our estimator for N{T), we use a Bayesian bootstrap technique 
that allows us to avoid the use of a specific — relation, getting rid of some possible 
bias thereby introduced. We therefore assume that the distribution of the estimator 
A^(> T) for a given density iV(> T) is of the following form: 



D{N)dN = d{\)d\ 



(6) 
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10"N(>4. keV) 

Fig. 6. The finite size of the sample implies that om estimate of the abundance of 
clusters has some uncertainty. The distribution of values that our estimator can take, 
given the sample we have, is estimated by a Bayesian bootstrap. Here we illustrate the 
histogram of the values taken by our estimator of the abundance of clusters at present 
epoch with temperature above 4 keV . A convenient analytical fit of this distribution is 
also shown: a ^ distribution with 6 degrees of freedom. 

with A = N/N{> T). We then use a Bayesian bootstrap to reconstruct the distribution 
function d: 10'* fake samples are built from the original sample, each having the same 
mean number of clusters, but with a given dispersion to account for Poisson noise. This 
procedure was used to find the distribution function of the estimator of iV(> T), as well 
as of the estimator for iV(> Ti) — A^(> T^)- The latter quantity is easier to handle with 
in a statistical analysis, as different temperatures can be chosen in such a way that the 
measurements are independent quantities (typically N{> Ti) — N{> T2), N{> T2)-'N{> 
),...). We found that a distribution function fits quite well the distribution function 
of the estimator, provided that the number of degrees of freedom is left as a free parameter 
(and is not necessary an integer) . An example is given in Figure |^. As one can see, the 
distribution provides an adequate fit to the overall distribution. From this we can 
infer confidence intervals on the density. The error bars given in Figure ^ reflect the 68% 
range of the value of the estimator of the density of clusters. 
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5. Cosmological Applications 
5.1. Mass-temperature relation 

In order to relate the observable properties of clusters to their mass, physical correlations 
are necessary. As discussed in the introduction, the temperature of the X-ray gas is 
probably the most convenient cluster mass estimator. Kaiser (1986) developed the scaling 
arguments which are essential to this approach (they were already used in 1980 in the 
pioneering by Pcrrenod, 1980). Some of the scaling arguments are exact, and they are 
remarkably powerful. In an Einstein de Sitter universe with initial fluctuations described 
by a power-law spectrum, the only process acting is gravity and therefore for the physics 
of the gas, scaling arguments will hold as long as no other mechanism plays an important 
role (like cooling or bulk heating). In such a situation, the only possible scale of the 
problem is associated with the non-linear scale, even in the presence of shock heating. 
This guarantees that the T — M relation for M, clusters follows an exact scaling law (for 
a given spectrum): 

nocM^/R*ocM^^^{l + z) 

It is reasonable to expect that such a relation will also hold between different masses 
forming at a given redshift and that is nearly independent of the power spectrum. How- 
ever, this remains an approximation as the geometrical aspect of the peaks corresponding 
to different masses at a given time may differ (d\ie to the difference in the height of the 
peaks). The same remark applies for two different spectra. Further checks are therefore 
needed. Analytical modeling of clusters adopting hydrostatic equilibrium is often com- 
pleted by the isothermal assumption. However, such a mass estimator can be misleading, 
even when temperature profiles are taken into accoimt (Balland & Blanchard, 1997). A 
further problem is the fact that a significant part of the pressure support of the gas could 
come from turbulent motions in the gas (Bryan and Norman, 1998b). Furthermore, the 
scaling laws are not expected to hold exactly in open models, significant departures may 
actually exist in this case (Voit & Donahue, 1998). It seems therefore safer to rely on 
numerical simulations in order to establish the mass-temperature relation. There are sev- 
eral ways to define the mass of clusters, which are not actually objects with well defined 
and sharp boimdaries. Considerable progress has been made in recent years thanks to 
Navarro, Frenk and White (1996), who showed that clusters are well fitted by a universal 
profile involving few parameters (which is not what would be expected based on the 
isothermal hydrostatic /3- model). Still, various definitions of the mass are used, result- 
ing in some confusion when one wishes to make detailed comparisons. The mass of a 
cluster could be defined as the mass within the so-called virial radius corresponding to 
a spherical region with a contrast density 6 = 178 (in the Einstein-de Sitter case). The 



Blanchard, Sadat, Bartlett & Le Dour: Cluster Temperature Function 17 

mass can also be defined hy 6 = 200, S = 500, or by the mass within a fixed physical 
radius such as the Abell radius or its equivalent in comoving coordinates. As most of 
the numerical checks of the mass function have been performed at the virial radius, it is 
certainly safer to estimate the mass at this very radius. The validity of the scaling models 
has been impressively demonstrated (Evrard 1989), and no noticeable difference has been 
detected among different spectra, masses and cosmological backgrounds. Furthermore, an 
intensive comparison of the results from various codes has verified that different numer- 
ical techniques do not lead to significant differences in inferred properties (Prenk et al., 
1999). This is especially true for the M — T relation. Even the inclusion of a significant 
energy injection does not introduce significant changes in these relations (Metzler and 
Evrard, 1998). The recent numerical work of Bryan and Norman (1998a) is certainly the 
most advanced in this area, and the scaling relations they obtained are impressive (see 
their figure 4) and convincingly independent of the underlying cosmological model. Their 
result can be summarized by the fact that the following mass-temperature relation: 

Ty = 3.8keyMi/^(OoA/178)^/^(l + z) 

(given for h — 0.5) reproduces well the results of their numerical simulations (M15 = 
M/10^^ solar masses). The statistical dispersion they found around this relation is only 
10%. The slope 2/3 might not be guaranteed to better than 15%. However, the normal- 
ization of the above relation in the various numerical simulations can differ (for instance, 
Evrard et al., 1996 found a value 20% higher than Bryan & Norman, 1998a, while the 
mean value obtained by Frcnk et al., 1999, is only 5% higher with a 5% dispersion among 
the various simulations) . The existence of turbulence on small scales (Bryan and Norman; 
1998b) indicates that one should be cautious on this issue, as the modeling of turbulence 
is often limited by the resolution of the simulation. Finally, it is important to notice that 
such masses arc 70% higher than what would be inferred from the isothermal hydrostatic 
equation (Roussel et al., 2000) and more than twice the mass inferred recently from de- 
creasing temperature profiles as considered by M98, or more recently by Nevalainen et 
al. (1999). In our modeling, we have adopted the above normalization, corrected for the 
effect of introducing a 10% dispersion in the M — T relation: 

Ty = 4.0keyMi/^(OoA/178)i/3(l + Z) 

{h = 0.5). We will examine the consequences of any change in this relation. 
5.2. Constraining the models 

The local temperature distribution function can be used to constrain the properties of the 
matter density fluctuations. Although this can also be applied to a non-Gaussian scenario 
(OB97; Robinson et al., 1999), we will restrict ourselves to Gaussian fluctuation models 
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Fig. 7. Contours giving the allowed range of values for the amplitude of matter density 
fluctuations on cluster scales. The continuous line is the 68% contour, the dotted line is 
the 90% contour, and the dashed line is the contour at 95%. 



only. This consists basically in finding the best parameters and confidence intervals from 
a fit to the temperature distribution function. However, there is a degeneracy between 
different cosmological models: equally satisfying fits can be obtained in either low- or 
high density models (OB92). Nevertheless, the density fluctuation amplitude can be 
usefully constrained for a given model. It is generally assumed that cluster abundances 
lead to a determination of the fluctuation amplitude at 8/i^^Mpc and therefore the 
normalization is communly expressed in term of erg. This is not completely correct. In an 
J7o = 1 universe, the temperature distribution function of clusters docs actually provide 
the normalization of the fluctuations on this scale. This comes from the fact that the mass 
of a 8/i^^Mpc sphere collapsing today will lead to a 4.5 kcV cluster, which is a typical 
temperature of a moderately rich cluster for which X-ray samples provide a reasonably 
well established abundance (at least this what the authors now believe!). In contrast, 
in a low-density universe (ilo ^ 0.3), a 8/i~'^Mpc sphere will produce a 2 keV cluster, 
for which we actually do not yet have reliable abundance estimates. This is not just a 
semantic difference, since for a low-density universe [flo ~ 0.2), the amplitude of erg 
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Fig. 8. Contours giving the range of allowed values for the index n of the primordial 
matter fluctuation spectrum on cluster scales for a Gaussian distribution function versus 
the normalization expressed in term of (jg (left side) and in term of Cc for a density 
parameter i7o = 0.4. The continuous lines are the 68% contours, the dotted lines are the 
90% contours, and the dashed hnes are the contours at 95%. 

could differ by up 50% (OB97). Following Blanchard & Bartlett (1998), we will use the 
following definition: 

and give constraints in terms of CTc. This has the additional advantage of being indepen- 
dent of the power spectrum: as erg and the power index n are correlated, the uncertainty 
on (78 is larger than on ctc- In order to illustrate this point we have plotted both quantities 
in Fig. ^. 

To find parameter constraints, we find the best model according to the following 
likelihood function: 

C = l[p{N{> T,,z,) - Ni> r,+i,z,)l^{f2o,n,-c}(> T.,,z,) - N[no..,aM> T.+i,z,)) (7) 

in which N is the measured abundance and -/V{no,Ti,(Tc} is value given by the model, 
which depends on the density parameter Qq^ the local slope of the power spectrum 
n and its amplitude specified by CTc- The probability distribution function p cannot be 
inferred from first principles and is therefore estimated with our bootstrap procedure. The 
maximum value of the likelihood has been normalized to unity, and confidence intervals 
on the parameters were estimated from the contours at £ = 0.6, 0.26, 0.13 corresponding 
respectively to the 68%, 90% and 95% contours (strictly speaking this holds only for 
normal distributions). The constraint on the shape of the spectrum n is rather poor. We 
find a broad range of possible values. The only obvious result is that models with fio = 1 
and n —1, representing the standard Cold Dark Matter spectrum, are excluded, as 
previously found (Blanchard & Silk, 1991; Bartlett & Silk, 1993). The effect of statistical 
uncertainties on the amplitude of matter fluctuations on cluster scales is quite small, and 
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therefore only a restricted range of values is allowed for a given Oq- The dependence on 
Oo is weak : 

(Jc « 0.6 1°;°^^ (95%) open case 

ac« 0.6 1°:°^^ (95%) flat case (8) 

(uncertainties are statistical and are given for fio = 1). Different values have been pub- 
lished for the amplitude obtained from the cluster normalization. These differences are 
mostly due to different cluster abundances used. Our normalization is slightly higher than 
what is often derived from abundances of X-ray clusters (but in better agreement with 
some normalizations based on optical surveys, e.g., Girardi et al, 1998). This is due to 
the fact that our number density of clusters at 4 keV is higher than previous estimations 
and that we used the Bryan and Norman (1998) normalization of the M — T relation 
which is slightly lower than that of Evrard et al. (1996). Finally, the differences with 
previous works concerning the dependence on fio come from the fact that we refer to CTc, 
which is the normalization on the true cluster scale rather than to as, which depends on 
the assumed spectrum (for n = — 1.5 we find erg = 0.96 in an open fio = 0.3 universe). 

It is interesting to try to evaluate the magnitude of possible systematic uncertainties 
on (Tc- As we have argued, the local cluster abundance now seems well established, so 
that we do not expect this to be a source of systematic error larger than the statistical 
uncertainty already quoted. The remaining dominant uncertainty comes from the use of 
the Press and Schechter mass function. Different results exist in the literature concern- 
ing its accuracy to describe the actual mass function (Lacey & Cole, 1994). Numerical 
simulations mostly suggest that the mass function is well described by the Press and 
Schechter formula: for instance, Borgani et al. (1998) have shown that it can be safely 
used for cluster modeling in arbitrary cosmology. 

5.3. Estimating flo 

The determination of the mean density of the universe is obviously a very important 
goal of observational Cosmology. Although it is fashionable to claim that current data 
favor a low value, it seems to us that the observational situation has not improved that 
much since the 1980's. Indeed, bulk flows on large scales lead to somewhat contradictory 
results and interestingly enough, fluctuations in the CMB seem to favor flat models 
and to exclude open ones (Lineweaver et al., 1997; Lineweaver and Barbosa, 1998; 
Le Dour et al., 2000). Perhaps the most stringent constraint comes from the baryon 
fraction observed in clusters (White et al, 1993). Still, a flrm upper limit on Oq is 
lacking, because of the uncertainties in primordial nucleosynthesis and cluster mass 
estimates (notice that virial masses estimated from the scaling models with the nor- 
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Fig. 9. Contours giving the range of allowed values for the index n of the primordial 
matter fluctuation spectrum on cluster scales for a Gaussian distribution function. The 
continuous line is the 68% contour, the dotted line is the 90% contour, and the dashed 
line is the contour at 95%. 

malization we adopt are nearly twice as large as those based on hydrostatic estimations!). 



The evolution of clusters is known to be a powerful probe of the mean density of the 
universe (Pcrrcnod, 1980; Peebles ct al. 1989). Evolution of the cluster abundance with 
redshift offers a new pass to flo (OB92, OB97) which is insensitive to the cosmological 
constant, and above everything else is really a global probe rather than a local one, as it is 
the case with the classical M/ L argument. It has been shown that the EMSS distribution 
can be modeled equally well within low- and high density universes, provided that some 
evolution in the Lx — Tx relation is adopted (OB97). The required evolution depends on 
fio as follows: 

L,aT3(l + ^)^ (9) 
with 

S ~ 4^0 - 3. (10) 

It is important to ask whether this modeling based on limited observational data 
is robust or not. Since this work, the observational situation has been considerably 
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improved, thanks to the ROSAT satelhte. The number counts for X-ray clusters are 
now available and fall well within the range of the predictions given by OB97. The 
ROSAT X-ray cluster luminosity function was derived from the BCS sample (Ebeling et 
al. 1997) and recently confirmed by De Grandi et al. (1999a). This observed luminosity 
function matches extremely well the one predicted by OB97. Finally, the amplitude of 
the fluctuations which was normalized to the original HA91 temperature distribution 
function, is in good agreement with what has been obtained in the present work. 
This considerably reinforces the strengt of OB97 modeling. Reichart at al. (1999) 
have recently re-analyzed the EMSS sample and have reached similar conclusions. 
Furthermore, similar value of B has been obtained by Borgani et al. (1998) from their 
analysis of the RDCS survey. The determination of Qq from the evolution of the — 
relation seems therefore possible, and we expect it to lead to reliable estimates. Thanks 
to the availability of ASCA temperature measurements for a significant number of 
high-redshift clusters, it has been shown that this relation undergoes at most only a 
mild evolution (in the positive direction, i.e., high-redshift clusters are more luminous 
than local ones with the same luminosity; Sadat et al., 1998). Using the relation between 
B and flo (eq. 9), Sadat et al. (1998) concluded in favor of a high value of the density 
parameter (CIq ^ 0.8) with a quoted uncertainty of 0.2 (this uncertainty only reflects the 
uncertainty due to the Poisson noise in the EMSS distribution, and does not take into 
account the uncertainty in the overall modeling). 

Thanks to the availability of the first X-ray selected sample of high-redshift clusters 
with measured temperatures (H97), it is now possible to constrain the density of the 
universe using directly the evolution of the X-ray cluster abundance to an intermediate 
redshift. A value of fio = 0.5 was obtained (H97) (this value has been abusively referred 
to as a low Hq.). Eke et al. (1998) who studied in detail the various systematic uncer- 
tainties have found a similar Qq- A higher value of f2o has been derived by Viana and 
Liddle (1999a) who discussed the importance of a proper correction from the uncertain- 
ties introduced by temperature measurement errors. Since these authors used a lower 
local abundance than the one we have inferred from the present work, it is interesting 
to re-address this question. Our primary goal here is to analyze the consequences of the 
higher cluster abundance at zero redshift that we derive. We leave to a future work a 
more complete investigation of the high-redshift sample. 

We have estimated the temperature distribution function of the distant cluster sample 
as given by H97, and estimated the suitable volume correction for the various cosmological 
models. Temperatures were corrected by the Bayesian term. The resulting temperature 
distribution function is given in Figure |o[ For comparison, we have also plotted the 
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Fig. 10. The integrated temperature distribution function at z = 0.33 inferred from 
Henry's sample is given by the continuous thick (red) Hne (corresponding to the case 
rio = 0.5). The two triangles are the estimation from Viana & Liddle (1999a, 1999b) 
with the original H97 sample. The line is the best fitting model. The continuous grey 
(green) line is the fit to our local temperature distribution function. The dashed grey 
(green) line is our best model at z = 0.33. 



temperature distribution function inferred from the local luminosity function. One may 
worry that our Bayesian correction is adequate. The bootstrap approach followed by 
Viana and Liddle (1999a) is certainly well adapted to the treatment of measurement 
errors. We have therefore compared their inferred abundances to ours and found very 
good agreement (using the 10 original clusters). This confirms that our analytic method 
provides a correction term which is as good as the bootstrap resampling technique. Our 
inferred temperature distribution is shown as the thick continuous line in Figure [l0| . 
We use the following likelihood function to infer the best-fit value of 

£ = []p (iV(> T,, Zi) - N{> T,+i,z,)\N{no,n,aM> ^r, z{) - N{n,.n,a^}{> T,+i,Zi)) 

i 

Y[p{Ni> T„Z,) - N{> r, + i,Z,)|iV{Oo,n,ae}(> r,,Z,) -iV{Oo,n,ae}(> Tj + l,Z,)) (11) 
3 

with = {3.,4.,5.4,6.6,8.,+oo,Zj = 0.05} and T, = {3.6, 4.4, 5.5, 6.7, 8., +cx), = 
0.33}. The distribution function of the measured abundances was estimated by a Bayesian 
bootstrap, as discussed in section 4.3. Because of the small number of clusters in the 
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Fig. 11. Figure a) represents the histogram of the values taken by our estimator from 
10'* bootstrap resampled samples (as in figure |6|). The spikes are real and reflect the fact 
that the abundance estimation is dominated by a single cluster. In this case, our analytic 
expression for the distribution function does not work very well (continuous line : 
distribution with 1.5 dof). b) represents the distribution function of the estimator with 
larger bins; in this case the fit is acceptable. 

sample, the inferred abundance looks very spiky (see Figure pT|a), which simply means 
that one cluster essentially dominates the statistic. In such a case, our fits are of course 
never as good as in the case where the number of clusters is larger, but they are reasonably 
acceptable when the binning in abundance is enlarged (see Figure ^T|b) . 

The likelihood function normalized to unity, is given in Figure |l^. The most likely 
values of fio ^^^e 0.92 (open case) and 0.865 (flat case). The shape of the likelihood 
function is very well fitted by a Gaussian, even if the probability functions we have used 
are significantly non-Gaussian. We can therefore make direct use of this function to give 
the confidence intervals at the la level. Our best fitting values are : 



no = 0.92i°J^^ (open case) 
r^o = 0.865l[JJ^5 (flat case) 



(12) 



This is very consistent with Sadat et al. (1998) value from a completely independent 
analysis. The constraint we have obtained is quite severe: an open model with CIq < 0.49 
is ruled out at the 95% confidence level, a fiat model with < 0.37 is ruled out at the 
same level, conclusions which are in clear disagreement with previous analyses based 
on the same high redshift sample (II97; Eke et al, 1998). The main differences come 
from our higher abundance at low redshift and from the fact that we explicitly take into 
account the effect of temperature measurement errors in the high-redshift sample. Note 
that Viana and Liddle (1999a) reached values which are consistent with ours. There are 
a number of other issues that differ in these various analyses and it is important to check 
whether these differences can result in significantly different values of JIq. Eke et al. 
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(1998) have analyzed various sources of systematics and concluded that they are not of 
critical importance. We reach similar conclusions for the effects they have investigated. 
For instance, we have modified the M — T relation by ± 20%, which changes Oq by ±0.1. 
The uncertainty in the validity of the Press and Schechter formalism is also significant 
but smaU. Using the revised version by Governato et al. (1998), we found values arround 
15% higher. We have checked that our likelihood approach is not biased by applying 
it to the predicted mean abundance in one model (fio = 0.85, n = —1.85, CTc = 0.606) 
and found that the best fitting model is (Oq = 0.85, n = —1.75, CTc = 0.61). We also 
have checked that varying the shape of the fitting expression does not change the final 
likelihood at any appreciable level. After this work was finished, we learned about 
Donahue & Voit (1999) and Henry (2000) works. The first authors found a significant 
evolution in the abundance of clusters over the redshift range [0.-0.8], in agreement 
with Blanchard & Bartlett (1998) and the present work. Nevertheless, they conclude it 
to be consistent with a lowdensity universe, but simultaneously they found a very fiat 
spectrum n ~ —2.2, which is excluded in our analysis based on the local temperature 
distribution function. Similarly, Henry (2000) found a moderate value of ~ 0.5 using 
the HA sample. Indeed, the strongest source of systematic error we found comes from 
the reference sample we used at zero redshift: using the abundances inferred from the 
sample limited at the bright end {fx > 4.10"^^ erg/s/cm^) leads to Qq = 0.65 consistent 
with Eke et al.(1998) and Henry (2000). With the fainter part of our sample (2.210"" 
erg/s/cm^ < fx < 4.10""'^^ erg/s/cm^) we obtained fig — 0.98. The difference obtained 
just by dividing the sample into two statistically equivalent sub-samples is surprisingly 
large (althouhg one may not worry of a difference at 1.5 a level): one would expect the 
uncertainty to come primarily from the high-redshift sample (which comprises only nine 
clusters) and not from the low-redshift sample. 

Our best fit value of ^lo inferred from the H97 sample is significantly higher than 
previous estimates from the same sample but it is consistent with the latest optical 
result by Borgani ct al. (1999). It is therefore important to examine the robustness 
of our analysis. This disagreement with previous results can be due to a higher local 
abundances of X-ray clusters, and a somewhat different treatment of the bias introduced 
by the errors at high redshift. As we have already argued, it seems unlikely that our 
local sample leads to a significant overestimation of the local temperature distribution 
function. In fact, one could argue that we may miderestimate the actual N{T), because 
of the possible incompleteness of our sample. However, regarding the agreement between 
our local temperature distribution function and what can be inferred from the luminosity 
function this is rather unlikely. The high redshift sample might be more worrisome. 
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As we have seen, a systematic difference may exist between two samples, leading to 
larger differences in flo than expected from Poisson noise. It is therefore conceivable 
that the high redshift sample is a statistical fluke. For instance. Eke et al. noticed that 
the temperature distribution within the redshift bin [0.3-0.4] in the original Henry's 
sample is statistically surprising as 9 of the clusters lie in the range [0.3-0.35] (now 8 are 
left). Another possibility is that high-redshift clusters are more massive than what one 
would infer from their apparent temperature (averaged over the luminosity which comes 
essentially from the core region). This could be possible for instance, if high-redshift 
clusters are more often dominated by cooling flows (or if the selection procedure favors 
cooling flow clusters), making them appear cooler than they actually are and producing 
an apparent evolution in the M — T relation. At the same time, in order not to produce 
an evolution in the L — T relation (not seen in the data) , these high-z clusters would have 
to be fainter. This is not very compelling, since cooling flows are expected to increase 
the X-ray luminosity. 

Another potential problem in the above determination of the density parameter is 
the quality of the EMSS sample: if the selection function is not well understood, it could 
be that the sample is missing a significantly larger fraction of the cluster population 
than expected. Indeed, clusters selection in EMSS is rather problematic, as the detection 
algorithm was designed to detect point sources. A mean correction for extended sources 
was applied (see Gioia and Luppino, 1994 and OB97 for more details), but one may 
nonetheless worry that this procedure is not well controlled. However, the deficit of 
high-z clusters observed in Figure ^ is of the order of 2 to 3. The possibility that the 
EMSS selection procedure could have missed clusters to such an extent seems unlikely. 
Furthermore, the modelling by OB97 predicted number counts which arc subsequently 
seen to be in good agreement with available ROSAT counts argues against significant 
incompleteness of the EMSS sample. An other potential problem could be a systematical 
bias in the EMSS fluxes. No evidence has been found by Nichol ct al. (1997). However, 
Ebeling et al. (1999) claim that a 40% offset in flux exists, which could explain half of 
the observed dimming. This possibility does not seem very appealing, because it would 
imply a very signiflcant evolution in the Lx — relation. In order to get rid of possible 
limitations of the EMSS sample, it will clearly be very important to see whether consistent 
results could be obtained from ROSAT selected samples of X-ray clusters. 

6. Conclusion 

The local temperature distribution function of X-ray clusters is an important tool for 
cosmology and can provide direct information on the dark matter distribution. It is 
therefore essential to have a good estimation of the local temperature distribution func- 
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Fig. 12. The likelihood function provided by the comparison of the abundances of clus- 
ters at z = 0.05 and z = 0.33. The continuous line corresponds to the open case, the 
dashed line corresponds to the flat case. The two horizontal line gives the la and 2a 
confidence interval on one parameter for a gaussian distribution. The dotted-dashed line 
is a fit by a two-tailled gaussian. 

tion. Furthermore, such a sample is of crucial importance as a reference sample in order 
to properly evaluate the evolutionary properties of the cluster population and for other 
cosmological applications. We have provided a new estimate of the temperature distribu- 
tion function for local clusters, based on a large sample of X-ray clusters with measured 
temperatures. This sample is essentially flux limited and we have argued that it is likely 
to be reasonably representative with a completeness estimated to 85%. We found an ap- 
preciably higher abundance of clusters than previous estimates, but in good agreement 
with the abundance inferred from the luminosity function and optical data. We have 
used the sample to study the statistical properties of the matter density fluctuations and 
obtained results consistent with previous works, although we obtained a normalization 
slightly higher, erg ^ 0.6 for ^Iq — 1. than previously found from X-ray clusters due to 
our higher abundance of 4 keV clusters. Probably the most important application of this 
new temperature function concerns the determination of the density parameter via the 
test of the evolution of the cluster abundance with redshift. In order to apply this test, 
we have used the Henry's sample which provides for the first time a direct estimation of 
the temperature distribution function at a non-zero redshift. We have found a clear indi- 
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cation that the abundance of clusters was smaller at the epoch corresponding to z ~ 0.33 
consistent with Donahue and Voit (1999). Using a likelihood approach, we inferred a high 
value of Hq ^ 0.9 ± 0.2. Low-density open universes {Hq < 0.49) are excluded at the 2 
a level. The exclusion region is nearly as severe in the flat case: fip < 0.37 models are 
excluded at the 2 a level. This result is entirely consistent with other independent analy- 
ses of the EMSS sample (Sadat et al., 1998; Reichart et al., 1999). We therefore confirm 
that the abundance of X-ray clusters as inferred from the EMSS favors a high-density 
universe. We have pointed out that this result is also consistent with what is known from 
existing ROSAT samples, although the situation is not as clear as in the case of the 
EMSS sample. It is important to keep in mind that the present method is one of the very 
few cosmological probes of ilo that is not based on local estimates, but is rather global 
in nature. However, given the importance of the conclusion, we believe that it should 
still be considered with caution. Our study of the local temperature distribution func- 
tion demonstrated that systematic uncertainties could be more important than expected. 
It is therefore essential to perform this test using an entirely different and independent 
sample. Temperature measurements with XMM of ROSAT selected clusters will allow to 
obtain such a sample. The application of the cosmological test will then probably lead 
to a more definitive conclusion concerning the mean density of the universe. 
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Table 1. X-ray temperatures, 90% uncertainties, ROSAT 0.1-2.4 keV luminosities (h= 
0.5) and fluxes. 

Ref. (1) Ebeling et al. 1996; (2) David et al., 1993; (3) Markevitch, 1998; (4) Fukazawa et 
al., 1998 (temperature are given with the emission from the central 0.05/i~^Mpc excluded); (5) 
Johnstone et al., 1998; (6) Arnaud & Evrard, 1999; (7) Arnaud, 1994; (8) Yamashita et al., 
1992; (9) Hughes et al., 1993; (10) Ebeling et al., 1998; (11) White et al., 1994; 



